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Abstract. The self-force acting on a (scalar or electric) charge held in place 
outside a massive body contains information about the body's composition, and 
can therefore be used as a probe of internal structure. We explore this theme by 
computing the (scalar or electromagnetic) self-force when the body is a spherical 
ball of perfect fluid in hydrostatic equilibrium, under the assumption that its rest- 
mass density and pressure are related by a polytropic equation of state. The body 
is strongly self-gravitating, and all computations are performed in exact general 
relativity. The dependence on internal structure is best revealed by expanding 
the self-force in powers of rj^^, with rg denoting the radial position of the charge 
outside the body. To the leading order, the self-force scales as r^'^ and depends 
only on the square of the charge and the body's mass; the leading self-force is 
universal. The dependence on internal structure is seen at the next order, , 
through a structure factor that depends on the equation of state. We compute 
this structure factor for relativistic polytropes, and show that for a fixed mass, it 
increases linearly with the body's radius in the case of the scalar self-force, and 
quadratically with the body's radius in the case of the electromagnetic self-force. 
In both cases we find that for a fixed mass and radius, the self-force is smaller if 
the body is more centrally dense, and larger if the mass density is more uniformly 
distributed. 



PACS numbers: 04.20.-q, 04.40.-b, 04.40. Nr., 41.20.Cv 

1. Introduction and summary 

An electric charge held in place in the curved spacetime of a massive body produces an 
electric field that responds to the spacetime curvature; a consequence of the interaction 
is a distortion of the field lines from an otherwise isotropic distribution near the 
charge, which leads to a net force acting on the particle. This is the physical origin 
of the electromagnetic self- force [T] , a subtle effect that is sensitive to the geometry of 
spacetime not just in the vicinity of the charge, but everywhere. 

Self-force effects in curved spacetime have been vigourously explored; for an 
extensive review, see Ref. [2]. Most of the recent activity was focused on the 
gravitational self-force, in an effort to model the inspiral and gravitational-wave 
emissions of a binary system with a small mass ratio [51 HI |S]. The prototypical 
problem, however, goes back to 1980, when Smith and Will [S] calculated the self- 
force acting on a particle with electric charge e held in place at a radius tq in the 
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Schwarzschild spacetime of a nonrotating black hole of mass M. The analogous 
problem of the self-force acting on a scalar charge q in the same spacetime was 
investigated by Wiseman ffl. The results are interesting: The radial compoment of 
the electromagnetic self-force (in Schwarzschild coordinates) is equal to e^Af/J^^/rQ, 
in which /o := 1 — 2Af/ro, while the self-force vanishes in the case of a scalar charge. 
The electromagnetic self-force is repulsive, meaning that the external force required 
to keep the particle in place is smaller when the particle is charged, compared to 
what it would be in the case of a neutral particle. These results are not intuitive: 
It is not obvious why the electromagnetic self-force should be repulsive, and it is not 
obvious why the scalar self-force should vanish; for an intriguing explanation based on 
Newton's third law, refer to Sec. IV of Ref. 0. 

The issue that interests us in this paper is the way in which the internal 
composition of the massive body affects the self-force acting on an electric or scalar 
charge. Because the self-force is the result of an interaction between the field and 
the spacetime curvature, and because this interaction will depend on the details of 
internal structure, it is expected that the self-force will reflect these details, and will 
therefore reveal some aspects of the body's composition. Our goal is to describe which 
aspects of the internal structure can be revealed by the self-force. 

This line of enquiry was initiated by Burko, Liu, and Soen [8^. Building on earlier 
work by Unruh [3], they calculated the self- force acting on electric and scalar charges 
held in place in the exterior spacetime of a massive thin shell. They found that their 
results agreed with the black-hole results to leading order in an expansion of the self- 
force in powers of r^^, but that structure-dependent terms appeared at higher-order. 
The cases of an electric charge outside a body of constant density, and of an electric 
charge outside a massive conductor, were examined by Shankar and Whiting (lOj . who 
confirmed the agreement with the black- hole results for large values of tq. 

The enquiry was recently pursued by Drivas and Gralla [llj . who considered 
static charges in the exterior spacetime of a massive body of arbitrary composition. 
They firmly established the universality of the self- force at order r^^ , and revealed 
a dependence on internal structure at order rQ^ and beyond. More precisely stated, 
Drivas and Gralla found that when the self- force is expanded in powers of r^^, the 
leading term that appears at order r^'^ is independent of the body's internal structure. 
In the case of an electric charge, it is always equal to the Smith- Will result e^M/rp, 
and in the case of a scalar charge, it is always equal to the vanishing Wiseman result. 
(This last statement holds for a scalar field that is minimally coupled to the spacetime 
curvature. For nonminimal coupling, the self-force continues to vanish when the charge 
is held outside a black hole, but it is equal to 2^q^M /rg when there is a material body, 
where ^ is a dimensionless coupling constant. In this case the universality of the scalar 
self-force at order r^^ partially breaks down: the self-force is only insensitive to the 
details of internal structure when the body is not a black hole.) 

The main concern of Drivas and Gralla was to display the universality of the self- 
force at order r^^ , and their other discovery, that there is a dependence on internal 
structure at order r^^ and beyond, was left largely unexplored. This is the theme 
that we intend to pursue here: How does the self-force depend on the body's internal 
structure, and which aspects of the internal composition can be inferred from a close 
examination of the self- force? These issues are interesting, because it is very difficult, 
in general relativity, to collect any information about a body's internal structure 
when the body is spherical, and when one is limited to external measurements. It 
is impossible, for example, to reveal any aspect of the internal structure by measuring 
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the motion of a test body. Going beyond test bodies, for example by measuring the 
orbital motion of a two-body system of comparable masses, reveals very little: the 
strong formulation of the principle of equivalence implies that the orbital motion is 
necessarily insensitive to the details of internal structure, until the bodies are so close 
together that significant tidal deformations are generated [12 HSl HH [13 [HI [H] ■ 

The self-force, on the other hand, can be used as a probe of internal structure. 
Indeed, one can easily imagine a thought experiment designed to exploit the self- 
force in this manner. Suppose that one is given a neutral particle, a second particle 
with an electric charge e (or a scalar charge q), and a massless rope to keep each 
particle in place in the body's gravitational field; the other end of the rope is held 
at a safe distance by an external agent. Keeping the neutral particle at a position rg 
requires the external agent to hold the rope firmly, and the tension in the rope reveals 
the body's gravitational force on the particle. Replacing the neutral particle by the 
charge produces a slight deficit in the tension, which corresponds to the repulsive 
action of the self-force; the self-force at position tq can thus be measured. Repeating 
the measurement for many positions allows one to infer the self-force as a function of 
ro, together with its expansion in powers of r^^; the terms of order Tq^ and beyond are 
those that reveal aspects of the body's internal stucture. There is, of course, nothing 
practical about this. But it is of interest to observe that as a matter of principle, one 
can infer elusive details of internal structure by careful measurements of a self-force 
acting on a static charge. 

To investigate the dependence of the self-force on the internal structure of a 
massive body, we consider a spherical star that consists of a perfect fluid in hydrostatic 
equilibrium. The fluid is characterized by a (baryonic) rest-mass density p and a 
pressure p, and these variables are taken to be related by the polytropic equation of 
state p = ifp^+^Z", in which X is a scaling constant, and n is the polytropic index 
(another constant). The star possesses a mass M and a radius R, and a (scalar or 
electric) charge is placed at a radius ro > -R in the vacuum region exterior to the star. 

Following Drivas and Gralla [TT], we compute the difference between two self- 
forces. The first is the actual self-force acting on the charge held in place outside the 
polytropic star, and the second is the self-force that would be acting if the charge were 
instead situated outside a Schwarzschild black hole; the bodies have the same mass, 
and the charge is placed at the same position outside each body. It is this difference 
that contains the dependence of the self-force on the body's internal structure. We 
show below that in the case of a scalar charge, the difference can be expressed as 

Af'j;alar--(|:) ^2(21 + I) Sr'""' Qiizo)Q',{zo) , (LI) 

in which zq := r^/AI — 1, Qi{zq) is a Legendre function of the second kind, a prime 
indicates differentiation with respect to zq, and Sf^^^^^ is a structure factor that 
depends on the stellar interior; an expression is provided below in Eq. (3.24). In 
the case of an electric charge we have instead 

3/2 



r, \ (zq - l)Q^(zo) 
£{i+l) 



Uzo), (1.2) 



,^0 

in which 5™ is the corresponding structure factor for the electromagnetic self-force; 
see Eq. (|4.21|) below. 
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Figure 1. Structure factor j^j. ^j^g scalar self- force, plotted as a function of 

R/M for selected polytropic equations of state labelled by the polytropic index n. 
The curve for a massive thin shell is also displayed for comparison. At large R/M , 
the structure factor increases linearly with R/M, with a slope that depends on the 
polytropic index. For some equations of state the structure factor is multi- valued 
at small R/M; the low-lying branch, however, corresponds to stellar configurations 
beyond the maximum mass, which are dynamically unstable. 



When ro ^ M the self-force differences are dominated by the leading term £ ~ 1 
in the sums, and the expressions reduce to 

af;,,,,,^ „5r'^'^^ (1.3) 

and 

4 e^M^ 

^F^r.--^Sr—^- (1.4) 

This is the statement that to leading order in an expansion of the self-force in powers of 
, the details of internal structure are revealed at order Tq^; and these are contained 
in the structure factors S'^^^i^'' and Sf^. 

In Fig. [l]we present plots of Sf^^^^'-' as a function of R/M for selected values of 
the polytropic index n. We find that for large values of R/M, the structure factor 
increases linearly with R/M, with a slope that depends on the polytropic index; these 
results imply that for fixed rg and M, the self-force increases linearly with the stellar 
radius. This is the same scaling that was found by Burko, Liu, and Soen j8j in the 
case of a massive thin shell, for which we also plot the structure factor. 

An increase of the self-force with the stellar radius is to be expected. The self- 
force results from an interation between the field produced by the scalar charge and 
the curvature of spacetime. Since the self-force would vanish in the pure Schwarzschild 
spacetime of a black hole, the interaction is limited to the region of spacetime occupied 
by the matter. For a fixed M, the region of interaction increases with the stellar radius 




Figure 2. Density profiles for selected polytropic models. The mass density 
p, normalized by the central density pc, is plotted as a function of r/R. Each 
model corresponds to a body of radius-to-mass ratio R/M ~ 15, and each 
curve is labelled by the polytropic index n. The corresponding values of 
the relativistic parameter b (defined in Sec. |2jl are 3.830543764152575e-2 for 
n = 0.5, 4.116715459636999e-2 for n = 1.0, 4.733363787126200e-2 for n = 1.5, 
5.977719194440032e-2 for n = 2.0, and 9.929978563520857e-2 for n = 2.5. The 
figure reveals that as n increases, the body becomes increasingly centrally dense. 



R, and it follows that the self-force should increase with R. This argument, however, 
does not explain why the enhancement with size is linear in R instead of some other 
relationship. The argument continues to apply to the case of a massive thin shell, for 
which the curvature is entirely concentrated within the shell; here also the size of the 
interaction region increases with the shell's radius, leading again to an increase of the 
self- force. 

In addition to the linear scaling of the self-force with the curves displayed 

in Fig. [l] reveal that for a fixed R/M, the self-force is a decreasing function of the 
polytropic index. An explanation for this effect can be deduced by extending the 
preceding argument to take into account the distribution of mass within the star. In 
Fig.[2]we show density profiles for polytropic models that share the same ratio R/M] 
one sees that models with larger values of n are more centrally dense than models 
with smaller values. A centrally-dense star will possess a smaller region of effective 
interation with the scalar field than a more uniform star, and this will produce a 
smaller self-force. 

In Fig. [3] we present plots of 5°™ as a function oi R/M for selected equations of 
state. We find that for large values of R/M, the structure factor increases quadratically 
with R/M , with a curvature that depends on the polytropic index; these results imply 
that for fixed tq and M, the self-force increases quadratically with the stellar radius. 
This is the same scaling that was found by Burko, Liu, and Soen in the case of a 
massive thin shell. 

The increase of the self-force difference with the stellar radius can be explained 
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R/M 

Figure 3. Structure factor 5j™ for the electromagnetic self-force, plotted as 
a function of R/M for selected polytropic equations of state labelled by the 
polytropic index n. The curve for a massive thin shell is also displayed for 
comparison. At large R/M, the structure factor increases quadratically with 
R/M, with a curvature that depends on the equations of state. 

in the same way as for the scalar case. Here the self-force difference results from 
an interaction between the electromagnetic field and the matter distribution, and a 
larger star gives rise to a larger region of interaction, and therefore a larger self-force 
difference; the argument does not explain why the enhancement with size is quadratic 
in R in the electromagnetic case. The property that the self-force difference, for a 
fixed R/M, is a decreasing function of the polytropic index, is also explained in the 
same way, with the help of Fig. [2j Here also the property is a consequence of the fact 
that stellar models with larger values of n are more centrally dense than models with 
smaller values. 

The main conclusion of this work is that aspects of the internal structure of a 
spherical star can be revealed in a close examination of the self-force at order r^^ . 
Our computations show that the self-force increases with R/M , in a linear manner for 
the scalar case, and in a quadratic manner for the electromagnetic case. In addition, 
we find that the rate of increase depends on the equation of state, a centrally dense 
body producing a smaller rate of increase than a more uniform body. 

In the remainder of the paper we establish the results presented previously. We 
begin in Sec. [5] with a description of relativistic polytropes. In Sec. [3] we describe the 
integration of the scalar-field equation in the spacetime of a massive body, and the 
computation of the self- force. In Sec. |4]we describe the calculations for the case of an 
electric charge. 



Self-force as probe of internal structure 



7 



2. Polytropic stellar models 

In this section we describe the polytropic steUar models that are involved in our 
examination of the influence of internal structure on the scalar and electromagnetic 
self-force. Our models are static and spherically symmetric, and we write the 
spacetime metric as 

ds^ = -e^'^' dt^ + f-^dr^ +r^{de^ +sm^ed(t)'^), (2.1) 

with V a function of r and / := 1 — 2m{r)/r. For a perfect fluid at rest the field 
equations are 

dm ,9 ,„ ^, 

— = 47rr^£7, (2.2) 

dr 

f =-(<'+>.)g = -^'(". + W.), (2.4) 

in which a is the energy density and p the pressure; the last equation is the condition 

of hydrostatic equilibrium. 

The polytropic equations of state are 

p = Xpi+i/", e = np, (2.5) 

in which K and n are constants, p is the rest-mass density, and e the internal 
thermodynamic energy; the total energy density is then a = p + e. 

To integrate the field equations wc introduce the scaling quantities pc (the central 
mass density), Pc '■= Kpl'^^^" (the central pressure), mo := 47rpc''o mass scale), and 
the squared length scale r^ := {n+ l)pc/(47rp^). We introduce also the dimensionless 
relativistic parameter 

b:=^= Kpl'^, (2.6) 

Pc 

in terms of which mo/ro = {n + l)pc/pc = (n + 1)6. We make use of the dimensionless 

variables ^, 9, and fi such that r = tq^, p = PcG^, and m = mo/x. In terms of the 
Lane-Emden variable 9 we also have p = Pc9"~^^, e = npc^""*"^, and a = pc{l + nb9)9^. 
The field equations become 

^ =e'(l+n60)r, (2.7) 
I = - [1 + (n + l)b9] (m + 6f r+i) , (2.8) 

with / = 1 — 2(n-|- l)6/x/^. The equations are integrated outward from ^ = 0, at which 
we impose the boundary conditions p, = Q, 6 = 1, and V = V'c- Integration stops at 
^ = ^1, at which p = pi, 9 = Q, and "0 = V'l- The star's radius is then R = ro^i, its 
mass is M = niopi , and the value of tpc is chosen so that e^^^^ = 1 — 2M/R, to ensure 
that the solution joins smoothly to the Schwarzschild metric when r > R. 

The numerical integration of the equations is facilitated by using v := p/^^ as a 
substitute for the mass function, and a; := In^ as a substitute for the radial variable. 
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Figure 4. Mass of a relativistic polytrope, in units of mo, as a function of R/M, 
for selected values of the polytropic index n. The lowest curve corresponds to 
n = 0.5, higher curves correspond to increasing values of n, and the highest curve 
corresponds to n = 2.5. In each case the mass reaches a maximum value at 
a minimum value of R/M. Configurations with smaller R/M are dynamically 
unstable and therefore unphysical. 



In terms of these we have 

^ = (1 + n66l)6i" - 3u, (2.10) 

ax 

rift P"^ 

^ = - y [1 + (n + 1)6^] {y + 6r+i) , (2.11) 

*^ = (n+l)6^(t/ + W"+i), (2.12) 
ax J 

with / = 1 — 2(n + 1)6^^1^ and ^ = e^. The boundary conditions are now placed at 
X = — oo, at which u = |(1 + nb), 9 = 1, and = ^c- In practice the integrations arc 
started at a; = a;o < such that = is very small, and starting values for y, 9, 
and V are obtained by expanding each quantity in powers of and determining the 
coefficients with the help of the differential equations. 

Integration of the field equations for a selected value of the polytropic index n 
gives rise to a continuous sequence of stellar models parameterized by b, which acts 
as a substitute for the central density pc- Because mo and ro depend on the central 
density (and therefore on b), it is necessary to rescale the mass and length units so as 
to eliminate this dependence before producing plots of the mass M and radius R as 
functions of b on the sequence. We therefore set 

M = mo6^^-"'/V, i? = fo6(^-"^/^a, (2.13) 

in which 
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are the rescaled mass and length units, respectively, which are independent of b and 
therefore constant on the sequence. Plots of M/toq as a function of R/M for selected 
values of the polytropic index n are presented in Fig. [4] 

3. Scalar self- force 

3.1. Scalar field and self-force in spherical spacetimes 

We consider the self-force acting on a static scalar charge q moving on a world line c in a 
curved spacetime; the world line is described by the parametric equations = z"(t), 
in which t is proper time. The charge produces a scalar potential $ that satisfies the 
wave equation 

^ -4T:q Jji{x,z{T)) dr, (3.1) 

in which □ :— .g^^VaV^ is the covariant wave operator, and S4{x,z) is a scalarized 
Dirac distribution. The self-force acting on the scalar charge is given by 

F°' =q{g"^ -\-u°'u'^)Vp^^, (3.2) 

in which u" :— dz°'/dT is the charge's velocity vector, and — $^ is the 

difference between $ and the Detweiler- Whiting singular potential [18] ; the regular 
potential is known to be smooth on c, and to be solely responsible for the self-force. 
For our particular application we take the spacetime to be static and spherically 



symmetric, we write the metric as in Eq. (2.1), and we choose c to be a static world 



line at r = To > -R. We decompose the potential in spherical harmonics, 

E 



$(r,0,0) = V$,„(r)r,^(0,0), (3.3) 



and substitution within the wave equation produces 



.d^^m^f^^r df di;\ d<^>eo ^(^+1) ^ . /2^+l 1/2^, . 

' ^+[^+2fdrr+''d^r^r —^^o = -47:q^^fo ^-'^o), (3.4) 

in which /o :— f{ro). Without loss of generality we placed the charge on the axis 9 = 0, 
and exploited the property Yfm(0, 0) = ^ {2£ -t- 1) /{4:n)Sm,Q of spherical- harmonic 
functions. The modes m 7^ of the scalar potential vanish. 

It is easy to show that the time and angular components of the self-force vanish in 
a static and spherically-symmetric situation, and that the radial component is given 

by 

= g/o lim ^[(9,$)^ - (9,$S)J ^ (3.5) 



in which 



i9r'^)i--- E ^ir)Yi^iO,^) (3.6) 



m— - 
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are the multipole coeffi cien ts of dr^, while {dr^^)i are those of the singular potential 



The limit in Eq. (3.5) can be taken by setting r ~ r^ + A, 6 ~ 0, and letting 



A — )■ (from either direction). With this choice, we find that 



The multiplole coefficients can therefore be obtained from the solution to Eq. (3.4). 
This equation must be integrated in the stellar interior of the polytropic model, and 
in the body's Schwarzschild exterior. 

3.2. Interior 

The equation satisfied by the mode of the scalar potential was given above in 



Eq. (3.4), and since the charge is situated in the exterior portion of the spacetime 
(ro > R), the right-hand side of the equation is zero everywhere within the interior. 
To integrate the equation we implement a change of variables from to rji defined 

by 

m ■■= —f, — ■ (3-8) 
The equation for becomes 

(-^ + m[vi-V + T^m J — = (3.9) 

when expressed in terms of rji and ^ r /ro, with 2? := 2 + ^f^^r{df/dr) + r{dip/dr) 
given explicitly by 



2?= 1 1^2- (n + l)6e'{2i^+ [l + (n- l)60]r}j 



(3.10) 



The equation is integrated outward from ^ = 0, at which ry^ = £. Here also the 
numerical treatment can be improved by adopting x = In^ as a new independent 
variable, and adding a few more terms to an expansion of rji in powers of 
Integration proceeds until ^ = ^i, at which rji = r]l. As we shall see, the number 
tjI is the entirety of the information required about the internal solution. 

3.3. Exterior 

Outside the matter the metric is described by the Schwarzschild solution, so that 



e^v = / = 1 - 2M/r. Equation (|3^ becomes 

+ 2(1 - ^).f - + 1,* . -*^,f^fi"Sir - r„,, (3,11, 

where we suppress the label "£0" on $. It is helpful to use z := r/M — 1 instead of r 
as a radial variable, and to re-express the differential equation as 

where zq :— rg/M — 1. Away from zq the right-hand side of the equation vanishes, 
and the solution is a linear superposition of Legendre functions Pi{z) and Qi{z). The 
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presence of the delta function implies that the solution must respect the junction 
conditions 



[$] = 0, 



dz 



= -47r- 



2£+l 



1 



(3.13) 



MM At: (zq - 1)i/2(^„ + i)3/2 

&i z = zq; here [/] := f{z — zo + e) — f{z — zq — e), with e approaching zero from 
above, is the jump of the function / across z = Zq. 



3.4. Matching 

The interior and exterior solutions must match smoothly at r = i?, or z = Z := 
R/M — 1, and the two exterior solutions (for z < zq and z > zq) must satisfy the 
junction conditions at z = zq. For z < Z we have that $ = $in, for Z < z < zq we 
have that $ = AP£{z) + BQp{z), where A and B are coefficients to be determined, 
and for z > zq we have that $ = CQi>{z), where C is a third unknown coefHcient; we 
exclude a term in Pe{z) when z > zq because it would produce a diverging field at 
z = 00. Letting a := APi{zo), (3 := BQi{zo), 7 := CQe{zo), and 6 := ^in{Z), we find 
that the matching conditions give rise to the relations 

a = (zo - l)(zo + l)Pi{zo)Qiizo) J, (3.14) 

P= - (zo - l)(zo + 1) g + ' %'%\ m^o)rJ, (3.15) 

[Z + l)Q'f {Z) -r]^Qi[Z)^ 

in which 

•^^=^"il/^(zo-l)i4„ + 1)3/2' (3-16) 

along with -f — a + /3 and S ~ [Pi{Z) / Pe{zQ)]a + [Qi{Z)/Qi{zo)]f3. To arrive at these 
results we used the Wronskian property of the Legendre functions, PiQ'^ — QgP'i = 
— (z^ — 1)~^, with a prime indicating differentiation with respect to z. As we shall see, 
the quantity of prime interest for a self-force computation is 7. 



3.5. Black-hole problem 

If the scalar charge were placed outside a black hole instead of in the exterior region 
of a spherical star, the Schwarzschild metric would apply everywhere instead of being 
restricted to r > i? or z > Z. In such a case there would be no interior solution to the 
scalar-field equation, and the two external solutions would be given by $ = APi{z) for 
z < Zq, and $ = CQi[z) for z > zq; a term in Qi{z) is excluded when z < zq because 
it would diverge on the event horizon. Letting a :— APi{zo) and 7 := CQi{zq), we 
find that the junction conditions at z = zq produce a — ^ = a. 



3.6. Self-force difference 

The radial component of the self-force acting on a scalar charge q at r = rp is given 



by Eq. (3.5), which we express as 



(3.17) 



with the understanding that the right-hand side is evaluated in the limit A — > 0. 
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In addition to this self-force, we may also consider the self-force acting on another 
scalar charge q placed at r = tq in the pure Schwarzschild spacetime of a black hole. 
For this situation we would have instead 



hole 9/0 [(^'■^holc)^ - 



(3.18) 



which refers to $hoic, the scalar potential as computed in the pure Schwarzschild 



spacetime, and to the same singular potential as in Eq. (3.17). The singular potential 
is the same in each case, because it depends on the structure of spacetime only in the 
immediate vicinity of the charge, where it is described by the Schwarzschild metric. 
Taking the difference between Eqs. ( |3.17[ ) and ( |3.18[ ), we have that 



hole ) 



(3.19) 



and we see that the difference between the two self-forces can be computed without the 
involvement of the singular field; this method of regularization was originally devised 
by Drivas and Gralla fTT^. Because the black-hole force actually vanishes [7], we have 
that 

FTtar = AF^ (3.20) 

and $hoie can in fact be identified with the Detweiler- Whiting singular potential. 
The force difference can now be expressed as 



2e+l 



dr 



-(ro + A) 



hole 



dr 



(ro + A) 



(3.21) 



in which each $ stands for the £0 mode of the spherical-harmonic expansion. While 
each individual mode-sum for F^^^^ and F^^^^ would fail to converge, the mode-sum 
for AF^ can be shown to converge exponentially. 

The mode-sum requires the evaluation of d^sta.T:/dr and d^i^oic/dr at r = rg -I- A. 
It is convenient to take A = 0+, so that r = tq is approached from above. We therefore 
insert Md^stai/dr = iQ'iizo) /Qi{zo) and Md<^\io\c/dr = lQ'i{zQ)/Qi{zo), and obtain 



2£+ 1 
47r 



(7-7) 



?£(2o)' 



(3.22) 



the sum now excludes £ = 0, because the monopole solutions are the same in the two 
spacetimes and therefore do not contribute to the difference. (The monopole solution 
describes the field of a spherical shell of charge g at r = tq; in both spacetimes the 
field vanishes inside the shell, and goes as g/r^ outside the shell.) Inserting the result 
obtained previously for 7 — 7 = /3, we arrive at the expression of Eq. ( 1.1 1, 



AFstar — — 





'z^~ 


t) 




\Zq + 





3/2 00 

^(2^-f i)5,g,(zo)g^(^o), 



in which 



{Z + \)P[{Z)-rf[P,{Z) 
{Z + l)Q',{Z)-r^jQe{Z) 



(3.23) 



(3.24) 



is a structure factor that depends on the stellar model (through the stellar radius R 
and the interior constant r/l) but is independent of tq. We recall that z := r/M — 1, 
Z :— R/M — 1, zq := ro/M — 1. and that a prime indicates differentiation with respect 
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Figure 5. Scalar self-force difference AF*", in units of (q/M)'^, plotted as 
a function of ro/M for selected polytropic equations of state labelled by the 
polytropic index n. Each stellar model corresponds to a body of radius-to- 
mass ratio R/M ~ 15 (see Fig[2]|. The fallo ff w ith tq is consistent with the 
asymptotic behaviour ~ r^^ predicted by Eq. (1.3 I, and for fixed ro/M the self- 
force difference is seen to decrease with increasing n; this behaviour is discussed 
and explained in Sec.[l] 



to z. The mode-sum can be evaluated straightforwardly once the numerical results 
for ril and R/M are available. 

The self-force of Eq. ( |3.23[ ) can be computed as a function of for selected stellar 
models characterized by a polytropic index n and a relativistic parameter b. Sample 
results are displayed in Fig. [5j and our findings were presented more fully in Sec. [T] 
The self-force for a massive thin shell can be computed by integrating the scalar-field 
equation in the flat interior and matching the solution to the Schw arzsch ild exterior; 
this yields rjj = Ij — 2M/R, which can be inserted within Eq. (3.241 and shown 
to reproduce the results of Burko, Liu, and Soen [8]. 



4. Electromagnetic self-force 

4-1. Electromagnetic field and self -force in spherical spacetimes 

We next consider the self-force acting on a static electric charge e moving on a world 
line c in a curved spacetime. The charge produces a vector potential that satisfies 
the wave equation 

n^a- r£<^,3 ^ -Aire ju°'5i{x,z)dT (4.1) 

in the Lorenz gauge V"$q = 0; here Rap is the spacetime's Ricci tensor. The potential 
gives rise to the electromagnetic field F^p — Vq$^ — V^^q, and the self- force acting 
on the electric charge is given by 

= ei^R^.w^ (4.2) 
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in which F^/^p :— — -F's'^ difference between the actual electromagnetic field 

and the Detweiler- Whiting singular field [THj; the regular field -Fr"^ is known to be 
smooth on c, and to be solely responsible for the self-force. 

For our particular application we take the spacetime to be static and spherically 



symmetric, we write the metric as in Eq. (2.1), and we choose c to be a static world 



line at r = ro > i?. The only relevant component of the vector potential is $t , and 
we decompose it as 

$t(r, e,4>)^Y. '^ttm{r)Yt^{e, cj,). (4.3) 

Substitution within the wave equation produces 

2Cp$tffl / r df dVA d^tm ^(^+1) ^ 
' +V+2fd-r- '^r^r 



= 47Te^'^^e^°f-'/^S{r-ro), (4.4) 

in which /o := /(fp) and Vo = 4'{i^o)- Once more we placed the charge on the 
axis 9 — 0, and exploited the property Yfm(0, 0) = \/ (2£ +1)/ {4:TT)Sm,o of spherical- 
harmonic functions. The modes m 7^ of the potential vanish. 



Inserting the decomposition of Eq. ( |4.3P within Eq. (4.2) yields 

- ee-^«/o lim - (9.$?)^, (4.5) 



in which 



i9r'^^\ - E ^W^-(M) = V^^(^o + A) (4.6) 

are the multipole coefficients of d^'^ti while {dr'^^)i are those of the singular potential 
$f . The limit a; — >■ z is taken by setting r — IS., 9 — Q, and letting A — > 0. The 



multiplole coefficients can be obtained from the solution to Eq. (4.4). As in the scalar 
case, this equation must be integrated in the stellar interior of the polytropic model, 
and in the body's Schwarzschild exterior. 

4-2. Interior 

The equation satisfied by the ^tm mode of the electromagnetic potential was given 
above in Eq. ( |4.4[ ), and since the charge is situated in the exterior portion of the 
spacetime (rg > R) , the right-hand side of the equation is zero everywhere within the 
interior. As in the scalar case we implement a change of variables from '^tm to 77^ 
defined by 

d\n<^>tm , . „x 

The equation for ^tto becomes 

^§+..(..-l)+I^..-'^=0, (4.8) 
d£. f 

with V := 2 + ^f^^r{df/dr) — r{dip/dr) given explicitly by 

V=j ^2- (n + l)6e'{4j^+ [l + (n + l)60]r}j . (4.9) 
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The equation is integrated outward from ^ = 0, at which r]£ ~ £. Here also the 
numerical treatment is improved by adopting x = In^ as new radial variable, and 
adding a few more terms to an expansion of r/t in powers of . Integration proceeds 
until ^ = ^1, at which ije = ijj. 

4-3. Exterior 

Outside the matter the metric is described by the Schwarzschild solution, so that 



g2V = J = 1 _ 2M/r. Equation (4.4) becomes 



)/^f^^'6ir^r,), (4.10) 



+ 2r/— - £{i + 1)$ = 4^(ee'^" 
ar^ dr 

I 1 /2 

where we suppress the label "t£0" on The factor e'^" is equal to /g , but we prefer 
to lump it with the electric charge e in order to keep the right-hand side in the same 



[$] = 0, 



ee^ 1 21+1 1 

mV 4^ (zo-l)i/2(zo + 1)3/2 



form as the right-hand side of Eq. (3.111; in this way the electromagnetic source is 
obtained directly from the scalar source by the replacement q —ee'^°, with ee'^° 
playing the role of an effective charge. 

Once again it is helpful to use z := r/M — 1 instead of r as a radial variable, and 
to re-express the differential equation as 

where zq :— ro/M — 1. Away from zq the right-hand side of the equation vanishes, 
and the solution is a linear superposition of (z — l)P^'(z) and {z — l)(5^(z), in a which 
a prime indicates differentiation with respect to z. The presence of the delta function 
implies that the solution must respect the junction conditions 

dz _ 

at z = Zq; we recall that [/] := f{z = zq + e) — f{z = zq — e), is the jump of the 
function / across z = zq. 

4.4- Matching 

The interior and exterior solutions must match smoothly at z — Z :— R/M — 1, and 
the two exterior solutions must satisfy the junction conditions at z = zq. For z < Z 
we have that $ = <i>in, for Z < z < 2:0 we have that $ = A{z- l)Pl,{z) + B{z - l)Q'g{z) , 
where A and B are coefficients to be determined, and for z > zq we have that $ = 
C{z — l)Q'f{z), where C is a third unknown coefficient. Letting a :— A{zo — l)P^'(zo), 
/3 := B{zo - l)Q^(zo), 7 := C{zo - l)Q^(zo), and S := $in(2'), we find that the 
matching conditions take the form of 

« = - ^'''^lj'j'^^^^^ nizo)Q'rXzo) J, (4.13) 

(zo-l)^(zo + l) eii+m{Z)-{l + 4){Z-l)P^{Z) 2 

^(£ + 1) ^(£+l)Q,(Z)-(l + r?i)(Z-l)Q^(Z)L^^^^°^J ^^■'^> 

in which 



ggVo 21+1 1 
■= -'"lirV An (zo-l)V2(zo + l)V2' (4-15) 

alongwith7 = a+/3and(5 = {[(Z-l)F;(Z)]/[(zo-l)P;(zo)]}a+{[(^-l)QK^)]/[(2o- 
l)(5^(zo)]}/3. As in the scalar case, the quantity of most direct interest is 7. 
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4-5. Black-hole problem 

If the electric charge were placed outside a black hole instead of in the vacuum region 
of a spherical star, the Schwarzschild metric would apply everywhere instead of being 
restricted to r > R oi z > Z. In such a case the two external solutions would be 
given by $ = A{z - l)Pe{z) for z < zg, and $ = C{z - l)Q^(z) for z > zq. Letting 
a :— A{zo — V)P[{zq) and 7 :— C(zo ~ l)Q£(-^o), we find that the junction conditions 
at z — Zq produce a = 7 = a. 



4-6. Self- force difference 

As in the scalar case we consider the difference between two self-forces, 



AF'' 



star 



hole 



ee 



/o 



hole) 



(4.16) 



the first acting on a charge e at position tq in the exterior of a material body, and the 
second acting on an identical charge at the same position outside a Schwarzschild 
black hole. Because the black-hole force is given by the Smith-Will expression 
-^hoio — e^AZ/p^^/ro [5, we have that the stellar force is given by 

--f^^' + AF\ 



star 



(4.17) 



Apart from a monopole piece that is entirely responsible for the black-hole force, $hoio 
can be identified with the Detweiler- Whiting singular potential. 
The force difference can be expressed as 



=0 



21+1 
47r 



d*s 



dr 



-(r-o 



A)-^(ro + A) 
dr 



(4.18) 



in which each $ stands for the £0 mode of the spherical-harmonic expansion of the 
vector potential $4. The mode-sum requires the evaluation of c?$star/c?^ and d$hoic/'^'' 
at r = ro -|- A. It is convenient to take A = O"*", so that r = tq is approached from 
above. We insert the appropriate expressions within AF*" and obtain 



ee 



fo 



M zo 



00 

e=i 



21+1 
An 



(7-7) 



'1 



£{£ + l)Qiizo) 
2^0-1 Q't(zo) 



(4.19) 



in which a prime indicates differentiation with respect to z :~ r/M— 1; the sum now 
excludes € = 0, because the monopole solutions are the same in the two spacetimes and 
therefore do not contribute to the difference. Inserting the result obtained previously 
for 7 



AF'' 



7 = /3, we arrive at the expression of Eq. ( 1.2 ) 

' Zq-I 



1 



3/2 00 

E 

£=1 



{2i + l)Se 



(zo-i)g^(zo) 
eii + 1) 



in which 



St 



i{i + i)Pe{z) - (1 + T^Diz - i)p;(z) 

£{e + l)QdZ) - (1 + f]}){Z - l)Q',{Z) 



Uzo), (4.20) 



(4.21) 



is a structure factor that depends on the stellar model (through the stellar radius R 
and the interior constant ry]) but is independent of tq. We recall that z := r/AI — 1, 
Z := R/M — 1, Zq := ro/Af — 1, and that a prime indicates differentiation with respect 




Figure 6. Electromagnetic self-force difference AF^, in units of (e/M)'^, plotted 
as a function of tq/M for selected polytropic equations of state labelled by the 
polytropic index n. Each stellar model corresponds to a body of radius-to- 
mass ratio R/AI ~ 15 (see Fig[2]|. The fallo ff w ith tq is consistent with the 
asymptotic behaviour ~ r^^ predicted by Eq. (1.4l, and for fixed ro/M the self- 
force difference is seen to decrease with increasing n; this behaviour is discussed 
and explained in Sec.[l] 



to z. The mode-sum can be evaluated straightforwardly once the numerical results 
for r]l and R/M are available. 

The self-force difference of Eq. (4.20) can be computed as a function of rg 
for selected stellar models characterized by a polytropic index n and a relativistic 
parameter b. Sample results are presented in Fig. [6j and our findings were discussed 
more fully in Sec. [T] As in the scalar case the self-force fo r a in assive thin shell can 
be computed by inserting 77] = i/y/l — 2M/R within Eq. (3.24); this reproduces the 
results of Burko, Liu, and Soen [8j. 
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